The Astrophysical Journal, 619: L127-L130, 2005 February 1 

Preprint typeset using I4TgX style emulateapj v. 6/22/04 
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ABSTRACT 

Polarization results from the Cosmic Anisotropy Polarization MAPper (CAPMAP) experiment are reported. 
These are based on 433 hr, after cuts, observing a 2 deg^ patch around the north celestial pole with four 90 
GHz correlation polarimeters coupled to optics defining 4' beams. The E-mode flat band-power anisotropy 
within I = 940;!;30Q is measured as 66;!;39 M^^i 95% confidence level upper limit for B-mode power within 
I = 1050!^f[! is measured as 38 fiK^. 

Subject headings: cosmology: cosmic microwave background — cosmology: observations 
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1. INTRODUCTION 

The cosmic microwave background (CMB) is arguably 
the most fruitful source of cosmological information. Its 
spectrum has been measured with extraordinary precision 
iFixsen et alJI 19971) and its spatial anisotropy has been very 
well characterized, with the Wilkinson Microwave Anisotropy 
Probe (WMAP; Bennett et al. 2003a) providing the most 
comprehensive results. There are two distinct patterns to the 
polarization of the CMB, conventionally termed E- and B- 
modes . The former have recently been dete cted at the degree 
scale ( Kovac et al.l l2002': Kog ut et alJl2003l) : they arise from 
the same density perturbations that dominate the temperature 
anisotropics. The latter, which can result from the lensing of 
E-modes due to the intervening matter distribution or from 
primordial gravity waves, are expected to be far smaller 
Characterization of the E-mode power spectrum is important 
for testing the understanding of the origin of the CMB as 
well as for breaking dege neracies in co smological parameter 
determinations (e.g.. .Bucher et alJ2001tiHu 2003^) . 

In this Letter we present results from the first season of 
CAPMAP (Cosmic Anisotropy Polarization MAPper), an 
effort to characterize the E-mode polarization at 90 GHz, 
where foregrounds are expected to be small, and at 4', where 
the polarization is expected to be near its maximum. In §§ 
|2]and|3]we describe the instrument, its characterization and 
calibration, and the observations. In subsequent sections we 
summarize the reduction of the data (§0, the analysis (§|5}, 
and the systematic studies (§|6j. We conclude this Letter with 
a discussion of the results and their implications (§Q. 
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2. INSTRUMENT 

The CAP MAP instrument is described in detail in 
iBarkats et al . (2005); here its most important features are 
given. In the winter of 2003, four W-band (84-100 
GHz) polarimeters (A-D) were in stalled in the foc al plane 
of the 7 m Cassegrain antenna (Ch u et al] Il978h at the 
Lucent Technologies facility in Crawford Hill, New Jersey 
(W74°ll'll", N40°23'31'). The horn and meniscus lens 
system attached to each polarimeter underilluminated the 
telescope mirrors with an edge taper of -35 dB on the primary, 
resulting in a final beam of 3'.88 FWHM (ctb = 1'.65) with an 
ellipticity less than 3%. The polarimeters were housed inside 
a single cryostat, nominally arranged so that A (C) observed 
0°.25 below (above) and B (D) 0°.25 left (right) of the optical 
axis of the telescope. 

Heterodyne analog correlation polarimeters that multiply 
two orthogonal components of the incident electric field 
produced an output proportional to one linear combination of 
the Stokes parameters Q and U. The two components derive 
from an ortho-mode transducer (OMT) oriented with one 
axis horizontal. The polarimeter output is thus proportional 
to the difference in power for field components aligned 
±45° to the horizontal in the focal plane. The first-stage 
low noise amplifiers (LNAs) in each arm are InP HEMT 
monolithic microwave-integrated circuits dev eloped by JPL 
and Northrop Grumman Space Technology JWeinreb et alJ 
1999). A mechanical refrigerator cools the horns, the OMT, 
and the LNAs to <40 K. One local oscillator inside the 
cryostat supplies the polarimeters with the 82 GHz to down- 
convert the 84-100 GHz signals to an intermediate frequency 
(IF) of 2-18 GHz. The oscillator signal in one arm of 
each polarimeter is phase-switched at 4 kHz, well above the 
1 // knee of the amplifiers. A detector diode monitors the 
total power in each arm. The phase-switched IF radiation is 
split into three 4 GHz- wide subbands (SO, SI, S2) prior to 
multiplication, yielding then 12 polarization and eight total 
power channels. 

The polarimeter outputs are digitized (24 bits, 100 kHz) 
with a SA analog-to-digital converter (ICS-610); digitally 
demodulated data are recorded at 100 Hz. The 100 kHz 
samples, averaged over 40 phase switches, are recorded to 
produce quadrature data sets, demodulated out of phase with 
the phase switch, which contain no celestial signal. 

3. OBSERVATIONS AND CALIBRATION 
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Between 2003 February 18 and April 6, 541 hr of CMB data 
were recorded. Calibration data were taken before, regularly 
throughout, and after the season. 

For the CMB observation the telescope's optical axis moved 
through the north celestial pole (NCP) with an 8 s period by 
±0°.53 at constant elevation. Polarimeter offsets from the 
optical axis are given in Table the sampled region is a 
circle of diameter 1°.6 centered on the NCP, with the inner 
circle of diameter 0°.14 missing. This scan strategy yields 
quite uniform Q/U coverage over the region. 




Fig. 1 . — The 20x72 data vector for receiver C. The scale ranges between 
±210 /^K; the average unchopped error for each element is 64 ^K. These 
pixels are spaced approximately 0.8 beamwidths in azimuth; in the LST 
direction, their spacing ranges from about 0.1 to 1 beam width depending 
on the azimuthal position. 



TABLE 1 
Polarimeter Position Offsets" 



A B C D 

co-elevation -0.022 -0.247 -0.025 0.200 
cross-elevation -0.286 -0.066 0.162 -0.059 

" The polarimeter offsets from the optical axis are given in degrees. For the 
CMB observations, the constant-elevation scan was centered with the optical 
axis at (-0°.03,-0°.09) from the NCP 

The pointing solution, determined with Jupiter 
observations, was refined to ±30" using 150 observations 
of astrophysical sources. The relative offsets, beam sizes 
and the total power calibration were determined from the 
observations of Jupiter in the total power channels. 

Noise temperatures were estimated from several elevation 
scans (20° - 90°). They ranged between 60 and 70 K and 
agreed with laboratory cold load tests. 

For the determination of polarized gains, a nutating plate 
was installed in front of the secondary mirror to produce a 
modulated polarized signal; this gave calibrations with an 
overall uncertainty of 10% and a 3% relative uncertainty 
among the channels. Frequent scans of the polarized 
source Tau A near parallactic angle 110° (where the signal 
in the receivers is maxim al) provided a good calibration 
check jBarkats e t al.l (20051). Polarimeter sensitivities are 
summarized in Table |2] 



TABLE 2 
Polarimeter Characteristics 





Polarimeter" 


<ff 


Channel 


A B C D 


SO 


3.2 3.4 2.9 3.1 


87 


SI 


2.1 2.1 2.1 2.2 


92 


S2 


2.7 2.0 3.5 2.5 


96 



" Sensitivities in mK^/s in thermodynamic units 
* Typical central frequency in GHz 

The calibration was corrected for slow few kelvin variations 
in the IF electronics and for a varying phase between the 
digitization and phase switch clocks. The IF gains changed by 
up to 3% K"' but the temperature was stable to 1 K for 70% 
of the data. The clock correction ranged from 0.9 to 1 and 
varied by at most 0.003 hr"'. The calibration for two of the 
channels (DSO, DS2) was adjusted by ^ 10% to compensate 
for gain compression. 

4. DATA REDUCTION 

The CMB data were divided into 40 s periods; data 
within each were divided into 20 ^ 3'-wide, nearly uniformly 



populated azimuthal bins, discarding the 9% of the data at the 
turnarounds. Means and standard errors calculated in each 
bin were carried forward in the analysis. Varying the period 
between 16 and 40 s did not alter the final results. 

A second-order polynomial was fitted to the azimuthal 
samples in each 40 s period and removed. The polarimetry 
channels reject the common mode intensity by more than 
23 dB, but their residual sensitivity to the atmospheric 
temperature, which varied daily by 10 K, dominated the 
intercept terms. The linear terms were evident when co-added 
over an hour, revealing slopes as large as 500 /iK deg"', with 
variations up to 200 /iK deg"'day"' . The quadratic terms were 
even smaller and more stable. 

Abnormally high fluctuations (> Aa) in the polarized 
channels or unusual changes in the ratios between different 
total power channels were observed for 20 hr of data. Both 
symptoms result from ice formation on the cryostat windows 
so these data were removed. The five most sensitive total 
power channels were used to identify bad weather with a cut 
on their fractional deviations. Only periods that passed the 
cut for all five channels were retained, eliminating 88 hr and 
leaving 433 hr for the final data sample. Ten percent variations 
in this cut level changed the result within statistics. 

Collapsing all polarization data into the 20 azimuthal bins 
reveals scan synchronous structures (SSSs) with rms levels 
from 10 to 26 /iK; from a 90 GHz side lobe beam map, these 
are likely due to ground pickup. Dependences on variables 
like system temperature, local sidereal time (LST), or ambient 
temperature were studied. Only BSl and BS2 showed such 
a dependence, varying with the ambient temperature. The 
mean of the data in each azimuthal bin was removed, and 
residual effects that could arise from correlations to ambient 
temperature are modeled as described below. The channel 
with the largest structure and variation (BS2) is dropped 
with negligible change in the results but with a significant 
reduction in the systematic uncertainty. 

After structure removal, corrections for atmospheric 
absorption (typically 20%) derived from the total power 
channels, were applied. The corrected periods, already 
parsed into azimuth bins, were binned into 72 LST bins, 
yielding 72x20=1440 element data vectors for each of the 
12 channels. The weighted averages of the three frequency 
channels (which have correlations less than 1%) were used in 
the likelihood analysis. One resulting data vector is displayed 
in Figure^ 

5. DATA ANALYSIS 

A Bayesian likelihood analysis was performed by 
maximizing the likelihood, 

/:=-=i=exp{-x^C-'x/2}, (1) 
Vdet(C) 
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where x is the data vector with 4 x 1440 entries and C is 
the covariance matrix. The vector x sums the signal {s) and 
receiver noise (n) terms: x, = AijSj + rii, where A,y determines 
the linear combination of Q (si) and U (sa) in pixel /. 
Similarly, the covariance matrix comprises signal (Ct) added 
to noise A^: C;y =< XjXj >= CTjj+Njj, where Nij = n]5ij, with 
typical noise n, ~ 60 /xK. The signal covariance depends on 
the multipole moments for E- and Z?-modes, Cf and Cf : 



(2) 



Here, the window function matrices account for (1) beam 
functions B] exp(-/^o'|), (2) the different pixel areas in the 
data vector, and (3) the conv ersion from (Q, U)to X = {E,B) 
iZaldarriaga & Seliaklll997i) . Figure |2] shows the trace of 
receiver C's ii-mode window function, similar to those of the 
other receivers. The likelihood is maximized with respect 
to one to three variables that parameterize Cf and/or Cf , as 
described below. 

The signal-to-n oise ratio eigenmode method was used (e.g., 
iBond et al.lll998h : transformations into the basis where the 
noise matrix is white and then into the eigenbasis of the 
transformed Ct are performed. Removing offsets every 40 s 
and subtracting SSSs transform x to Dx. Accordingly, C is 
transformed to DCD^, and the analysis is continued in the 
subspace onto which D projects. 

The code was exercised extensively with simulated data 
sets. Realizations with noise levels from that of the data to 
20 times smaller were superposed on CMB realizations of 
the concordance model ( S pergel e t al. 2003). The likelihood 
estimator was found to be unbiased; however, at the noise 
level of the data, the loss of sensitivity from the offset 
removals is a factor of 2. 

Four likelihood evaluations were made: first by taking Cf 
as a multiplier of the concordance model, with Cf = 0; second 
by taking 1(1 + l)Cf to be a constant (flat band-power), with 
0; third by simultaneously taking Cf as a multiplier to 



the concordance model and l{l+l)Cf as a constant; and fourth 
by taking Cf to comprise three flat band-power levels, in the 
/ bins (2-500/501-1500/1501-2900), with Cf = (three-band 
analysis). Fit results are s ummarized in Tab le Effective Is 
are calculated according to lSteinhardtl (ll 995h and quoted with 
the central 68% region of Wi /I. 

The various likelihood evaluations yield a 2 ct detection 
of the E-mode signal, with no B-mode signal or significant 
dependence on frequency observed. Results from the three- 
band analysis are shown in Figure |2 The upper band has 
negligible correlation with the other two. For the lower 
two bands the marginalized results from integrating over 
their likelihood contours are displayed. The marginalization 
shifted the peak of the middle band down by 6%, increasing 
its width by 12%. A fit of the middle /-band's likelihood to an 
offset lognormal distribution ( Bond et al..200 0) gives 66 /iK- 
with a variance of (44 /iK^)^ and a noise related offset of 
93 /iK^. 

6. EVALUATION OF SYSTEMATIC UNCERTAINTIES 

The three most important sources of systematic uncertainty 
were associated with the telescope pointing, the scan 
synchronous residual structure and the relative gains. Their 
impact was evaluated on fits of E-mode power to a multiplier 
of the concordance model; selected studies confirmed that 



CAPMAP r 
Window \ 





Fig. 2. — Likelihood results. Top: Results from the three-band analysis, the 
concordance model (long dashed line), and the trace of the window function 
, as defined in equation |2|, for receiver C, in arbitrary units (dotted 
line). For the outer bands, 95% confidence level upper limits are displayed. 
Lower left: Individual marginalized likelihood curves for the three bands. 
Lower right: Concordance model fitted to the data and several null maps as 
described in the text. 



uncertainties estimated with those fits apply to the three other 
fitting methods. 

The effect of possible mispointing was determined by 
varying the global pointing solution within its uncertainty, 
~ 10% of a beamwidth. The 1 ct shifts in azimuth and 
elevation were added in quadrature, resulting in an uncertainty 
of 12%. 

Effects of the SSSs were simulated. It was determined 
separately in each channel for two halves of the data, divided 
according to ambient temperature and its extracted (linear) 
variation. Multiple simulations were made with structure 
added to the CMB according to the ambient temperature in 
each period. This shifted the result upward by on average 0. 1 
(times concordance), and the spread, as the underlying CMB 
was varied, was ±0. 1 . The uncertainty is then conservatively 
taken to be ±0.2 (corresponding to a 10% uncertainty). 

The effect of relative gain uncertainties was estimated 
from the scatter in the result with the gain of each of the 
11 polarimeters chosen randomly within the estimated ±3% 
uncertainty, leading to an uncertainty of 6%. 

Three other sources of uncertainty were studied. The 
beam size was uncertain to 2%; varying it by even 10% 
leads to less than a 3% change in the result. As mentioned 
earlier, the polarimeters have a weak response to temperature 
anisotropics (-23 dB), and the optics produces a small 
quadrupole response at scales smaller than the beam (-12 dB). 
Simulations show that these also lead to negligible effects (see 
alsolHu et alJ2003h. 

Several checks of the robustness of the result were 
performed. Excluding single radiometers or 25% data 
samples gave consistent results. Difference maps were 
created by splitting the data into two parts according to high 
and low system and ambient temperatures and time. Data 
vectors were produced for each half (Di and D2) and the 
likelihood was evaluated on (Di-D2)/2. These likelihoods, 
together with that for the quadrature data, are shown in 
Figure |2l none show a signal. Finally, a parallel analysis 
was performed. Offsets were removed on a 4 s timescale, 
and an independent code implemented a different pixelization 
scheme and a different means of evaluating the associated 
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covariance matrix and nulling out the appropriate modes. It 
was used to obtain a multiplier to concordance and gave a 
consistent result. 

Adding the three main uncertainties in quadrature, the 
overall systematic uncertainty becomes 17%. The 10% 
overall gain uncertainty contributes an additional 20% 
without reducing the significance of the detection. 



TABLE 3 
CAPMAP RESULTS" 

ii-mode multiplier to concordance, 1=1120^^^^, Cf =0 



frequency 
result* 



all" 
2.0-^ 



SO 

3 9+'-3 



SI 

-0 1+2** 



S2 

Q 4+8.8 



Fits to flat band-power, Cf =0 except for fit for B 



1 

result" 



1050!«« 
< 38AiK^ 



1120!«o 
31!39^k2 



El" 
410+™ 



E2" 
940+330 

^ -300 



E3" 
1940^°° 



" Simultaneous fit with B also gives 2.0+j j 

* With 68% interval of highest posterior density or 95% upper limit 
^ Fitted simultaneously with E 
Three-band analysis 



7. DISCUSSION 

Evidence of polarization anisotropy in the CMB in the 
multipole region around / = 1000, where the E-mode power 
spectrum is expected to peak, has been presented. Signals 
of increasing significance are seen when fitting for (1) a flat 
band-power, (2) a multiplier of the expected concordance 
power spectrum, and (3) power in a band centered at these 
high Z -values. The E-mode signal does not change when 
fitting simultaneously for both E- and B-modes, and the latter 
are not detected. Here the possible contaminations to these 
results are discussed together with their significance. 

Foregrounds are not expected to be significant. From 
Ifinkbeiner et al. (1999) dust maps and WMAP synchrotron 
maps (iBennett et al.l l2003bl) around the NCP, temperature 



anisotropies of 4.4 and 1.3 fjK, respectively, are expected. 
Even with pessimistic assumptions about polarization 
fractions and how these foregrounds extrapolate to 
CAPMAP's angular scales, their contributions are well 
below the observed signal. 

There are no deep source maps exac tly at the NCP. A 
survey at 90 GHz (Holdawav et alJll994 extrapolates 4400 
sources brighter than 100 mJy over the entire sky, implying 
a 20% chance for one to be in the CAPMAP field. From 
the instrument response (3 fiK mJy~'), such a source would 
contribute at the 20 /iK level in 1 pixel if it had the same 
high polarization as Tau A (7.5%). Effects on the polarization 
power spectra, along with those from weaker but more 
copious sources, are negligib le. Conclusions f rom a dedicated 
source survey near the NCP (Leitch et al. 2000) are similar. 

Adding the systematic uncertainty in quadrature to the 
result in the middle 1-band, a null result is excluded at 
just over the 2 a level. Including the overall calibration 
uncertainty, this result is consistent with the concordance 
model within 72% of the highest posterior density. This 
result at 90 GHz is consistent with the very recent 
results from the Degree Angular Scale Interferometer and 
Cosmic Background Imager at 30 GHz JLeitch et al.ll2004t 
iReadhea cjet a l. 2004) . 
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